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Abstract 

Deconfinement and chiral symmetry restoration are explored in a confining, 
renormalisable, Dyson-Schwinger equation model of two-flavour QCD. An or- 
der parameter for deconfinement is introduced and used to establish that, in 
the chiral limit, deconfinement and chiral symmetry restoration are coinci- 
dent at T c ~ 150 MeV. The transitions are second order and each has the 
same critical exponent: (3 ~ 0.3. The deconfinement transition is found to 
exhibit sensitivity to the current-quark mass. f n and m v change by no more 
than 10% for T < 0.7 T c , however, as T — > T c , thermal fluctuations cause the 
pion bound state contribution to the four-point quark-antiquark correlation 
function to disappear. 
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The Dyson-Schwinger equations [DSEs] provide a nonperturbative, renormalisable, con- 
tinuum framework for analysing quantum field theories. An important example is the 
fermion DSE, which has proven useful in the study of confinement and dynamical chiral 
symmetry breaking [DCSB] [|TJ. The DSEs form a tower of coupled equations, which must be 
truncated to arrive at a tractable problem. Truncations that preserve the global symmetries 
of a field theory are easy to implement. Preservation of a gauge symmetry is more difficult 
but there is progress in that direction H. The elementary quantities are the Schwinger 
functions whose analytic properties provide information about confinement and DCSB; e.g., 
the absence of a Lehmann representation for the two-point dressed-gluon Schwinger function 
(gluon propagator) entails the absence of asymptotic gluon states; i.e., gluon confinement. 
The approach is reviewed in Ref. |§. 

The finite temperature properties of QCD are important in astrophysics and cosmology, 
and may be explored in a future-generation heavy-ion accelerator programme. Theoret- 
ical tools that can be employed reliably in the nonperturbative study of deconfinement 
and chiral symmetry restoration are therefore valuable. In Ref. 0] a one-parameter model 
dressed-gluon propagator was proposed, which provided a good description of n and p-meson 
observables. The single parameter in this model is a mass scale, mt, that marks the point 
where the nonperturbative, infrared enhancement found in gluon DSE studies || becomes 
dominant. The model gluon propagator has no Lehmann representation. The calculated 
quark propagator also has no Lehmann representation; therefore both the gluon and quark 
are confined. With m t fixed at T = 0, one has a renormalisable DSE-model of QCD, which 
manifests both confinement and DCSB, whose finite temperature behaviour may provide 
insight into the finite temperature properties of QCD. 

The DSE for the renormalised dressed-quark propagator at finite temperature involves 
a sum over Matsubara frequencies. In nonperturbative DSE studies, where the analytic 
structure of the dressed-quark propagator is calculated rather than assumed, it is necessary 
to perform this sum numerically because the usual analytic methods of evaluating it rely 
upon the quark propagator having a Lehmann representation, which may not be the case 
for a confined quark. The absence of a Lehmann representation also complicates, if not 
precludes, a real-time formulation of the finite temperature theory. 

Herein we employ a Euclidean space formulation with {7^,7^} = 25 nu , 7 M = 7^, j3 = 74, 
7-p = J2i=i liViy u n = (2n + 1)ttT and fl n = 2niiT, where T is the temperature . The DSE 
for the renormalised quark propagator, S = [— • p&a(p, u n ) — i/3 u n ac(p, oj n ) + <tb(p, <^n)]> 
is 



S (p, Ul n ) 




(1) 



mbm is the bare mass and the regularised self energy is 

27 • p E' A (p, u n ) + i/3 u n E' c (p, u n ) + Tf B (p, uj n ) , 



(2) 



with 




(3) 
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where T = _A,B,C; V A = -(Zf /p 2 )i 7 • p, V B = Z u V c = -(Zjujip, and j£ g = 
T J2m=-oo f A (C)3- m Eq. (§|), T u (q,uj rn ;p J uj n ) is the renormalised dressed quark-gluon 
vertex and D tiU (p,Q n ) is the renormalised dressed-gluon propagator. One can also write 

S~ l (p,uj n ) = 27 ■pA(p,u n ) +i(3u n C(p,uj n ) +B(p,u n ). (4) 

In renormalising we require that 



p 2 +ui 2 =/j, 2 



ij ■ p + i(3 uj + m R 



(5) 



which entails that the renormalisation constants are given by: Z£(p, A) = 1 — ujq] A), 

Z 2 (/i, A) = 1 — T,' c (p,uq;A), mn(fi) = Z 2 mb m (A 2 ) + T,' B (p,uo;A), and the renormalised self 
energies are 



T{p,u) n \[i) = K T + T,' :F (p,uj n ;A) - E^(p, w ; A) , 



(6) 



JF = A, B, C, = 1 = Kc an d Kb = mji(p). 

It is invalid to neglect A(p, w n ; p) and C(p,u> n ;p) and their dependence on their argu- 
ments. In studying confinement the (p, a;„)-dependence of A and C is qualitatively important 
since it can conspire with that of B to eliminate free-particle poles in the quark propaga- 
tor 0]. Further, all calculated observables are quantitatively sensitive to A and C; e.g., in 
the present study we find A(0,l} ) ~ 1.5 and C(0, luq) ~ 1.5 — 1.8 (increasing with T) and 
they return to 1 at the renormalisation point - neglecting this variation entails f n — > 2f n , 
which one may not be able to compensate by readjusting a model's parameters. 

To study the chiral limit one must solve the DSE with w,r{p) = 0, not = 0. We 
denote solution functions obtained in the chiral limit by a subscript "0"; e.g., B (p,u n ; p) 
and a Bo (p,u n ;p). 

The finite-T gluon propagator is 



g'D^p, ft) = P L (p, Q)A F (p, ft) + Pi(p)A G (p, ft) 



in Landau gauge, where 



Plip) 



0; /i and/or v — 4, 



ft 



p2 



1,2,3 
= 1,.. 



(7) 



(8) 



,4. A "Debye-mass" for the 



gluon appears as a T-dependent contribution to A F . 

We employ a parameter-free finite-T extension of the model dressed-gluon propagator of 
Ref. @: A F (p, ft) = V(p, ft; m D ) and A G (p, ft) = £>(p, ft; 0); 



X>(p, ft; m 

Y 



m 2 S 0n S 3 (p) - 



-(p 2 +n 2 +m 2 )/(4m 2 ; 



ft 2 



(9) 



where d = 12/(33-2JV» and the "Debye-mass" is ml = cT 2 , c = 4vr 2 dc, c = (N c /3 + N f /6), 



which is included in a manner analogous to that in Ref. 
in A G . 



]. There is no T-dependent mass 
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The first term in Eq. (^ is an integrable, infrared singularity || that generates long- 
range effects associated with confinement ||. The second term ensures that, neglecting 
quantitatively unimportant logarithmic corrections, the propagator has the correct pertur- 
bative behaviour at large spacelike arguments. In this model the large-distance, confining 
effects associated with 5o n 8 3 (p) are completely cancelled at small distances and one recovers 
the perturbative result. T>(p, fl; m) has no Lehmann representation and hence represents a 
confined particle, since this ensures the absence of gluon production thresholds in 5-matrix 
elements describing colour-singlet to singlet transitions. 

Equation (^) is an Ansatz based on T = studies augmented by finite-T perturbation 
theory. Solving the gluon DSE at finite-T would provide additional insight but no such stud- 
ies are currently available. Without further information; for example, solving the coupled 
gluon-quark DSE system or lattice estimates of the two-point gluon Schwinger function, we 
cannot be sure that the T-dependence we have introduced provides a qualitatively accu- 
rate representation of the temperature evolution of the gluon propagator. The paucity of 
relevant experimental data means that one cannot presently use experiment to constrain 
this extrapolation. This is the primary source of uncertainty in our study and we therefore 
advocate a cautious interpretation of our results near the transition temperature. 

The model of Ref. || is recovered with two further simplifying specifications. The ap- 
proximation 

Z x = Z 2 and = Z^ (10) 

is used. Gauge invariance in finite-T QED entails these Ward identities. In QCD the 
analogous identities are more complicated and this is only a simplifying truncation, which 
has proven qualitatively and quantitatively reliable in Landau gauge, T = studies. In 
nonperturbative DSE studies p|,|l0| Z 2 and Z 2 are finite and ~ 1. 



In addition, the rainbow approximation 

r M (g,u; m ;p,u; n ) = 7^ (11) 

is used. In T = studies this has proven to be reliable in Landau gauge; i.e., an effica- 
cious phenomenology with a more sophisticated Ansatz only requires a small quantitative 
modification of the parameters that characterise the small-A; 2 behaviour of the gluon prop- 



agator [H,|ll 



Confinement can be investigated by studying the analytic properties of the two-point 
Schwinger function of a given excitation. The confinement test proposed in Ref. ||, 
and used to good effect in Ref. is appropriate to this task. Consider Ab (x, 0) = 
(T/27rx)££lo AS^s), 

2 r°° 

a b () ( x ) = ~ d PP sin(px) a Bo (p, u) n ) . (12) 

71 JO 

For a free fermion of mass m: &b{p, w n ) = m /[ LjJ n + P 2 + rn 2 ], hence A B (x) = 



m exp(— x yto 2 + m 2 ). This illustrates that the n = term dominates the sum. 

For a free fermion, M(x; T) = —(d/dx) In A^ o (x) = ^(tt 2 T 2 + m 2 ), a finite constant. 
Finite-T effects become important for T ~ m/ir. In DSE studies, m is a mass-scale charac- 
teristic of dynamical mass generation: M^j d m 300 MeV, hence we expect finite-T effects to 
become noticeable at T ~ 100 MeV. 
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A Schwinger function with complex-conjugate poles, with §l(p 2 ) < (timelike) and 
Q{p 2 ) oc b, does not have a Lehmann representation; i.e., it represents a confined excitation 
with no associated asymptotic state. This is the nature of the quark Schwinger function 
obtained in Ref. [||. In such cases A° Bq (x) has at least one zero. Denote the position of the 
first zero in A° Bo (x) by n, and define uq = 1/ri, then deconfinement is observed if, for some 
T = T c , Uo(T c ) = 0; at this point thermal fluctuations have overwhelmed the confinement 
scale-parameter and the poles have migrated to the real-axis. 

We solve the DSE for S(p,u n ) numerically with 



mt = 0.69 GeV , m R (ji) = 1.1 MeV 



(13) 



which were fixed in Ref. by requiring a best x 2 -fit t° a range of 7r-meson observables at 
T = 0. We use d — 4/9 in Eq. @ and renormalise at \i = 9.47 GeV ||. With these choices 
our current study has no free parameters. We use A/fi = 1, |^ nmax |/|p| max ~ 1, n\ p \ = 64 
points in the |p|-array. Our results are cutoff independent. At T = 5 MeV the results in 
Table I of Ref. HJ] are reproduced to within 6%. 
At finite-T the pion mass is given by 



mlNl 



(m R (fi) (qq)f,)n 



(14) 



with (/ dux 



(2T)/(vr 2 ) E^ = o Jo A d PP 2 
(m R (n) (qq)ti)n = 
AN r f dfi T B (a Bo - 



Bn 



which vanishes linearly with m R (fi), and 

Nl = N c fd^ T Bl{ 



o \ 



2 2 i 2 2, 2 



'C) 



(15) 



+ °B(r' B ] ~\v 2 ([<£ [oc<y" c ~ 
+ P 2 (oao"a - (°'a) 2 ) + °bv" b - ^'b) 2 ])} 



(16) 



with a' B = d<JB{p 2 ,uj n )/dp 2 etc., is the canonical normalisation constant for the Bethe- 
Salpeter amplitude in ladder approximation. The pion decay constant is obtained from 



f N 

J 7V 1 * 7V 



2N C J dfi T B [a A a B + \p 2 (a' A a B - cr A a' B )} 



(17) 



Equation (|14[) provides an accurate estimate of the mass obtained in solving the dressed- 
ladder Bethe-Salpeter equation for the pion Equations (|15D -( p!7|) neglect corrections of 
order m 2 , which are unimportant (< 3%) for all values of T; i.e., m n is dominated by a 
linear response to m R at all T. 

In deriving these formulae one assumes = ij 5 B , where is the pion Bethe-Salpeter 
amplitude, which is a manifestation of Goldstone's theorem in the DSEs Since B is 
nonzero only if chiral symmetry is dynamically broken, these formulae are not valid above 
any chiral symmetry restoration temperature, T c x . 
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We employ the simplest order parameter for DCSB: 



X = B (p = 0,u ). (18) 

We plot X (T) and u {T) in Fig. § The curves, fitted on 120 < T (MeV) < 150, are of 
the form a (1 — T/Tc) 13 with 

a x = l.lGeV, /3 X = 0.33, 2? = 150 MeV . . 

a U0 = 0.16 GeV , (3 Uo = 0.30 , T»° = 150 MeV . ' 

The transitions are coincident and have the same critical exponent, within errors: ~ 10%. 

The massive quark is not deconfined at the same temperature. However, the behaviour of 
v m undergoes a qualitative change at T"° = T x , where no contribution to B(p,ou n ) remains 
that does not vanish as m R -> 0. For T < T x : a Vm = 0.15 GeV, (3 Vm = 0.23, T" m = 
160 MeV; for T > T x : a Um = 0.26 GeV, (3 Um = 0.71, T£ m = 180 MeV. For T > T x , a fit 
with j3„ m = 0.23 has a five-times larger standard deviation. The u/d-qusak in this model 
is therefore deconfined at T ~ 180 MeV, suggesting that the deconfinement temperature 
increases with current-quark mass. 

We find that, on 120 < T(MeV) < 150, N%, N n / w and (m iJ (/i)(gg) M ) are described by 
a (1 — T/T C )P, with T c = 150 MeV in each case and 

a N 2 = (0.18 GeV) 2 , (3 N 2 = 1.1, 
a UN l = (0.15 GeV) 2 , f3 h Z = 0.93, (20) 

a (m R {qq)) = (0.15 GeV) 4 , (3(m R (qq)) = 0.92. 

In order to avoid round-off errors associated with division near T = T c we used these results 
to determine 

a m7r = 0.12 GeV, /3 Wvr = -0.11, 
a u = 0.12 GeV, (3 fn = 0.36, 



with T™- = 150 MeV = T/-. We plot m^(T) and f n (T) in Fig. @. 

We explored the finite-T properties of a renormalisable, confining, DSE-model of 
QCD Q. Introducing an order parameter for confinement we found that, in the chiral limit, 
a deconfinement transition at T « 150 MeV is accompanied by the coincident restoration of 
chiral symmetry. The single model parameter, fixed at T = 0, is appropriate for two-flavour 
QCD: the transitions are second order and are not described by mean-field critical expo- 
nents. Similar results have been obtained in recent numerical simulations of lattice-QCD at 
finite-T [0. 

/,!■ and m-n are weakly sensitive to T for T < 0.7 T c x . However, as T approaches T c x , the 
mass eigen value in the pion Bethe-Salpeter equation moves to increasingly larger values, 
as thermal fluctuations overwhelm attraction in the channel, until at T = T c x there is no 
solution and f n — > 0. This means that the pion-pole contribution to the four-point, quark- 
antiquark correlation function disappears; i.e., there is no quark-antiquark pseudoscalar 
bound state for T > T x . This may have important consequences for a wide range of 



physical observables [14], if borne out by improved studies; e.g., this T-dependence of f n 



and m n would lead to a 20% reduction in the it — > jiv^ decay widths at T « 0.9 T x . 
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FIGURES 




T (MeV) 

FIG. 1. The parameters for the curves, x{T) (diamonds) and vq{T) (circles), are presented in 
Eq. (H). 
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FIG. 2. f w (T) (diamonds) and m 7r (T) (circles) plotted against T. See the discussion around 
Eq. (P). 
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